Pauli Tomography: complete characterization of a single qubit device 



F. De Martini*, Giacomo Mauro D'Arianot, Andrea Mazzei* & Marco Riccf 
*Istituto Nazionale di Fisica della Materia, Dipartimento di Fisica 

Universita "La Sapienza" , Roma, 00185 Italy, 
' Istituto Nazionale per la Fisica della Materia and Dipartimento 
di Fisica "Alessandro Volta" 
Universita di Pavia, Via Bassi, 6, 1-21100 Pavia, Italy 



The marriage of Quantum Physics and Information Tech- 
nology, originally motivated by the need for miniaturization, 
has recently opened the way to the realization of radically new 
information-processing devices, with the possibility of guaran- 
teed secure cryptographic communications, and tremendous 
speedups of some complex computational tasks. Among the 
many problems posed by the new information technology [I.L. 
Chuang & MA. Nielsen, Quantum Information and Quantum 
Computation, (Cambridge Univ. Press, Cambridge, 2000).] 
there is the need of characterizing the new quantum devices, 
making a complete identification and characterization of their 
functioning. As we will see, quantum mechanics provides us 
with a powerful tool to achieve the task easily and efficiently: 
this tools is the so called quantum entanglement, the basis of 
the quantum parallelism of the future computers. We present 
here the first full experimental quantum characterization of 
a single-qubit device. The new method, we may refer to as 
"quantum radiography", uses a Pauli Quantum Tomography 
at the output of the device, and needs only a single entangled 
state at the input, which works on the test channel as all pos- 
sible input states in quantum parallel. The method can be 
easily extended to any n-qubits device. 

How do we usually characterize the operation of a de- 
vice? Actually, we are interested just in linear devices, 
since quantum dynamics is intrinsecally linear. Any lin- 
ear device, either quantum or classical (examples are: an 
optical lens or a good amplifier), can be completely de- 
scribed by a transfer matrix which gives the output vector 
by matrix-multiplying the input vector. In quantum me- 
chanics the inputs are density operators pi n and the role 
of the transfer matrix is played by the so called quantum 
operation M of the device, that here we will denote by 
E. Thus the output state p ou t is given by the quantum 
operation applied to the input state as follows: 



Pout 



E (pin) 

Tr [E{p m )] 



(1) 



and the normalization constant Tr[E(pi n )] is also the 
probability of occurence of the transformation E, when 
there are other possible alternatives, such as when we 
consider the state transformation due to a measuring de- 
vice for a given outcome. 

Now the problem is: how to reconstruct the form of 
E experimentally? One would be tempted to adopt the 
conventional method JlJ] of running a basis of all possi- 
ble inputs, and measuring the corresponding outputs by 
quantum tomography^. However, since the states p arc 
actually operators, not vectors, in order to get all possible 



matrix elements we would need to run a complete orthog- 
onal basis of quantum states \n) along with their linear 
combination 2 a (|n') + i k \n")), with k=0,l,2,3 and i de- 
noting the imaginary unit (this is a simple consequence 
of the polarization identity). However, the availability of 
such a set of states in the laboratory is, by itself, a very 
hard technological problem (states with a precise varying 
number of photons and, even worst, their superposition, 
are still a dream for experimentalists). 

The quantum parallelism intrinsic of entanglement 
now comes to help us, running all possible input states in 
parallel by using only a single entangled state as the in- 
put! This was first shown in 0. Hence, we don't need to 
prepare a complete set of states, but just one entangled 
state, a state commonly available in modern quantum 
optical laboratories! 



I. PAULI TOMOGRAPHY 

Assume for semplicity, and with no loss of generality, 
that the entangled state spans a 4-dimcnsional Hilbcrt 
space Hi ® H2, each space H^ of dimension d=2. The 
input entangled state, 



(2) 



expressed in terms of the basis vectors |n) (£> \m) = \nm) 
of the two spaces H^ (i =1,2) is used as shown in Fig. 1 
where the entangled systems consist of two single-mode 
optical beams. The key feature of the method implies 
that only one of the two systems, say system i=l, is in- 
put into the unknown device E, whereas the other is left 
untouched. This setup leads to the output state R otl t, 
which in tensor notation writes as follows 



R out = E® / I 



(3) 



where I denotes the identical operation. It is a result of 
linear algebra that R «t is in one-to-one correspondence 
with the quantum operation E, as long as the state | 
is full-rank, i. e. it has non-vanishing components on 
the whole state-space of each system, such as, for in- 
stance, a so called maximally entangled state. Full-rank 
entangled states can be easily generated by Spontaneous 
Parametric Down Conversion of the vacuum state, as in 
the experiment reported here. Note that by this method 
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the problem of availability of all possible input states is 
solved: we just need a single entangled state | Vt')), which 
works as all possible inputs in a sort of quantum paral- 
lelism! 

Now, how to characterize the entangled state R on t at 
the output? We obviously need to perform many mea- 
surements on an ensemble of equally prepared quantum 
systems, since, due to the no-cloning theorem [|| we can- 
not determine the state of a single system ||. For this 
purpose a technique for the full determination of the 
quantum state has been introduced and developed since 
1994. The method named Quantum Tomography |j| has 
been initially introduced for the state of a single-mode 
of radiation, the so called Homodyne Tomography, and 
thereafter it has been generalized to any quantum sys- 
tem. The basis of the method is just performing mea- 
surements of a suitably complete set of observables called 
quorum. For our needs, we just have to measure jointly 
a quorum of observables, here the spin observables cr^ 
(i=0,l,2,3), on the two entangled systems at the output, 
in order to determine the output state R ou t, and hence 
the quantum operation E. 

In this paper we present the first complete experimen- 
tal characterization of a quantum device, which in our 
case will be a single-qubit device. The qubit is encoded 
on polarization of single photons in the following way 



\0) = \l) h \0) v 
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namely with the "logical zero" state corresponding to 
a single horizontally polarized photon and the "logical 
one" state corresponding to a single vertically polarized 
photon. In the following we will denote by h and v the 
annihilation operators of the horizontally and vertically 
polarized modes of radiation associated to a fixed wave- 
vector, k. Using single photon states we encode a qubit 
on the polarization. In the polarization representation, 
the Pauli matrices write as follows: 



o-i 



/iV 



o~2 = i (Wv — v^ti) 
0-3 



■ 1 T 1 



(5) 



The ring of Pauli matrices is completed by including the 
identity do = h'h+v'v. In the following we will denote by 
a the column three- vector of operators a = (ui, 01, 03), 
and by a the column tetra- vector a = (<To, o~i, 02, 03), and 
use Greek indices for three- vectors components a =1, 2, 
3 (or a = x, y, z), and Latin indices for tetra- vector 
components: i = 0, 1, 2, 3. 

A wave-plate changes the two radiation modes accord- 
ing to the matrix transformation: 



Wi 



(6) 



where the matrix W(0, 9) is given by: 
WU,6) ■■ 



cos a — sin f 
sin 9 cos 9 



1 

e** 



cos sin t 
— sin 9 cos ( 



z + + cz_ sz_ 

SZ- Z + — CZ- 



(7) 



where s = sin 29, c = cos 29, 9 is the wave-plate orienta- 
tion angle around the wave- vector k, z± = \ (l ± e l< ^) , 
cf) = A is the wave-lenght and 5 is the optical path 
through the plate. Special cases are the j plate which 
can be used with 9 = ? to give the right and leftg circu- 
larly polarized modes 



' ' X =w(- - , , 

1 2' 4) V v 



/4 ( h + iv 



y/2 \-ih + v 



(8) 



and the -| plate wich can be used to give the diagonal 
linearly polarized modes 



W 



J_ ( h + v 

V2\ h - V 



(9) 



The Heisenberg picture evolution of the Pauli matrices 
corresponding to the unitary U on the qubit Hilbert space 
will be given by a rotation, which we will denote as fol- 
lows 



E/W = R(U)a. 



(10) 



In particular, for a 0— wave-plate we have the rotation 
matrix 

s 2 +c 2 cos(j) — c cos (f) sc( 1 — cos 4 
R(W((j),9))= I csin0 cos0 — ssin</> 
sc(l — cos0) ssin0 c 2 + s 2 cos< 



(11) 



In particular, for a -|-wave-plate we have 



W(ir,9) = 



cos 29 sin 29 
sin 29 -cos26> 



- cos 46> sin 49 
R(n,6) = I -1 

sin 49 cos 49 



(12) 



The <7 z -photo-detector is achieved as in Fig 2. From Eq. 
pi] we can see that a u^-detector can be obtained by pre- 
ceding the a z -detector with a A/2-wave-plate oriented 
at 9= tt/8, whereas a er^-detector is obtained by preced- 
ing the <7 z -detector with a A/4- wave-plate oriented at 9= 
tt/4. When collecting data at a <7 a -detector, we will de- 
note by er Q = ±1 the random outcome, with a a = ±1 
corresponding to the h-detector flashing, and a a = — 1 
corresponding to the v-detector flashing instead. The ex- 
perimental averages for the complete setup must coincide 
with the following theoretical expectation values 



- (1) 4 2) = ((^KC/t®/)^ 1 ' ®af ] ){U® J) I*)) (13) 



(i) 



s i^ s< o^ an< l 



(1) (2) 
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and, in particular, s 

(n) 

s^ denoting the random outcome of the detector of the 
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nth beam (n=l,2) in the entangled state. For maximally we obtain the matrix ^ nm for the input states in terms 



-,(!) 



si 2) = for all 



entantangled states we have also s a 
a = x,y, z. The theoretical expectations for the setup 
without the device to be characterized are be given by: 

((* | <8> af) |*» = Tr [vP+a^a*] (14) 

Where Vp denotes the matrix of the state on the custom- 
ary basis of eigenvectors of a z . In particular, for the four 
Bell states 



1 



(15) 



we have 
1 



V2 



H = 



II- 11 

III- 1 
1-1-1 1 



1-11 



1 



(16) 



Using Eq. [fi] and || we obtain 

^~s% = J2 R ^(U)&Xf,(V) 



(17) 



In particular, in the lab we use the multiplet state corre- 
sponding to $ = <r x /-\/2, whence we have 



44 = R a p(u)^(^) 



II. METHOD, APPARATUS AND RESULTS 

The input maximally entangled state |\P)), expressed 
by Eq 2, was SPDC generated in the laboratory by an op- 
tical parametric amplifier (OPA) physically consis-ting of 
a nonlinear (NL) BBO (fi-barium-borate) crystal plate, 2 
mm thick, cut for Type II phase matching and excited by 
a pulsed mode-locked ultraviolet laser UV having pulse 
duration r = 140 fsec and wavelength (wl) A p =397.5 nm. 
Precisely, the apparatus was set to generate on the two 
modes k^, i.e. the entangled systems i =1,2, single pho- 
ton couples in a polarization entangled "triplet" state, 
viz. |*)) = 2" 1 / 2 \a x )), according to Eq. 15. The wl 
of the emitted photons was A= 795 nm. The measure- 
ment apparatus consisted of two equal polarizing beam 
splitters PBSi with output modes coupled to four equal 
Si-avalanche photo-detectors SPCM-AQR14 with quan- 
tum efficiencies QE ~ 0.42. The beams exciting the de- 
tectors were filtered by equal interference filters within 
a bandwidth AA = 6 nm. The detector output signals 
were finally analyzed by a computer. We want now to 
determine experimentally by this apparatus the matrix 
elements of the state | VP)) expressed by Equation 2. This 
can be achieved as follows: from the trivial identity 



(nm\ VP)) = VP r , 



of the following ensemble averages 



VP, 



((vP| 01) (nmj vP)) 

Tpj oi) (oi| W) 



(19) 



where the unmeasurable phase factor is given by: 
exp(i<p) — Vtoi/ |vPoi| • The choice of the vector 1 01) is 
arbitrary as it is needed only for the sake of normaliza- 
tion: e.g. we could have used 1 10) or |11), instead. Using 
the tomographic expansion over the four Pauli matrices 
H, we see that, in virtue of Eq. |l^, the matrix el- 
ement of the input state is obtained from the following 
experimental averages: 



4^5 



Qij(nm)sf\f' 



(20) 



where 



P=«*|01)(01| *)) = i(l + 4 1 )) (l-4 2) ) (21) 

is the fraction of coincidences with both cr z -detectors fir- 
ing on h, and the matrix Q(nm) is given by 



/ (nm) 



10) 



ID 



(22) 



the unitary matrix U nrn of the device is now obtained 
with the same averaging above but now for the state at 
the output of the device: | f/VP)) = (U <g) I) | VP)) . There- 
fore, we now have: 



(U*) r 



((f/\P| 01) (nm| vP£/)) 



(23) 



where we use again Eqns.^l|, |2^, but now the average ex- 
pressed by 20 is carried out over the output state | C/VP)). 
The (complex) parameters U ram are obtained from Eq. 
p3| by matrix inversion. This is of course possible since 
the matrix ^ is invertible, in virtue of the entangled char- 
acter of | VP)). 

The experimental demonstration of the tomographic 
process is given in Figures 3 and 4 where both real and 
imaginary parts of the four components of the matrix U 
are reported for two different " unknown" devices inserted 
in the mode k! of the tomographic apparatus. The exper- 
imental results are shown together with the correspond- 
ing data evaluated theoretically. Furthermore the exper- 
imental "variance" of the data are also reported. Each 
"unknown" device is represented in one figure, namely: 

Figure 3: a single Waveplate [<p = 0.457r; 9 = — 0.1387r] 
i.e., with retardation phase tp — (0.457r) and orientation 
angle respect to the "horizontal" direction "h" : 9 = 
(-0.1387I-). 

Figure 4: a combination of 2 Waveplates: a Waveplate 
[(p = 0.45tt; 9 = -0.138tt] followed by a A/2 Waveplate 
[tp = tt;8 = +0.29tt]. 

As we may see, the experimental results are found in 



(18) good agreement with theory. 
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III. CONCLUSIONS 



Figure Captions 



We have given the first demonstration, in a simple 
singlc-qubit context, of a novel Tomographic method 
which is able to fully characterize the properties of any 
device acting on a quantum system by exploiting for the 
first time the complete intrinsic parallelism of the quan- 
tum entanglement. This method establishes a new fun- 
damental framework of utterly paradigmatic relevance in 
the domains of modern Quantum Measurement theory 
and Quantum Information. 

Our method is expected to be of general and far reach- 
ing relevance. In facts, it can be adopted within more 
general and complex multi-qubit systems. For instance 
by this method a full characterization of a two-qubits de- 
vice, such as a controlled-NOT, can be achieved. In this 
case we just need to double the input and the measure- 
ment setup, by providing two input entangled states and 
four detectors coupled at the outputs of the two pairs 
of the device channels. The full quantum characteriza- 
tion of the the device is finally obtained by a joint tomo- 
graphic reconstruction on both channels of the device. 

We are indebted with the FET European Network 
on Quantum Information and Communication (Con- 
tract IST-2000-29681-ATESIT), the I.N.F.M. PRA 2001 
"clon" and with M.U.R.S.T. for funding. 
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FIG. 1. General experimental scheme of the method for 
the tomographic estimation of the quantum operation of 
a single qubit device. Two identical quantum systems, 
e.g. two optical beams as in the present experiment, are 
prepared in an entangled state | V I')). One of the systems 
undergoes the quantum operation E, whereas the other is 
left untouched. At the output one makes a quantum to- 
mographic estimation, by measuring jointly two observa- 
bles from a quorum {0(1)}. In the present experiment 
the quorum is represented by the set of Pauli operators. 

FIG. 2. Pauli-matrix measurement apparatus for photon 
polarization qubits inserted at the end of each test optical 
beam. 

FIG. 3. Experimental characterization by Pauli Tomog- 
raphy of a single 
optical Waveplate [</? = 0A5ir;8 — — 0.1387r] inserted on 
channel ki with the following optical properties: retarda- 
tion phase if = (0.457r); orientation angle of the optical 
axis respect to the laboratory horizontal direction "h": 
6 = (— 0.13871"). The experimental real and imaginary 
parts of the four matrix elements Uy of the Waveplate 
are shown together with the related measured statisti- 
cal variances. The corresponding theoretical values are 
shown for comparison. 

FIG. 4. Experimental Chractcrization by Pauli Tomog- 
raphy of a combination of two optical Waveplates: the 
Waveplate [tp = 0.45tt; 9 = -0.138tt] (cfr.Fig.3) followed 
by a A/2 Waveplate [<p = tt;9 = +0.297r]. The experi- 
mental real and imaginary parts of the four matrix ele- 
ments Uij of the combination are shown together with 
the related measured statistical variances. The corre- 
sponding theoretical values are shown for comparison 
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